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ABSTRACT 

An optical vortex, produced at one point in an optical beam, would propagate through an optical system to 
another point where the vortex can be used for some purpose. However, asymmetrical optical elements in such 
a system can cause astigmatism or at least distroy the rotational symmetry of the beam, which may affect the 
propagation of the vortex in an undesirable way. While an optical vortex in a rotationally symmetric, stigmatic 
Gaussian beam retains its initial morphology for as far as it propagates, the morphology of an optical vortex 
in an asymmetric or astigmatic Gaussian beam changes. The vortex can even be replaced by another with the 
opposite topological charge. We consider the behavior of single noncanonical vortices propagating in Gaussian 
beams that are asymmetric and/or astigmatic. General expressions for the vortex trajectories are provided. The 
locations of the flip planes and the evolution of the anisotropy of the vortex are considered for different non-ideal 
situations. 
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1. INTRODUCTION 

One of the significant advances that is being made in nano-technology is the ability to manipulate small objects 
and particles, even to the level of single atoms. A highly focused optical beam can be used as an 'optical tweezer' 
with which single atoms can be moved around. On an even more sophisticated level, optical beams can also be 
used to rotate or twist small particles. In the early 90's Allen and coworkers discovered that an optical beam 
can carry orbital angular momentum 1 (in addition to spin angular mometum, which is associated with circular 
polarization). This angular momentum can be transferred to a small particle that absorbs the light. 

Optical orbital angular momentum is associated with phase singularities, called optical vortices, which were 
discovered by Nye and Berry in the 60's. 2 An optical vortex is a point on the cross section of a beam (or a 
curved line moving along the beam). At this point the intensity of the beam is zero. Around this point the phase 
increases by an integer multiple of 2-7T. This integer is called the topological charge of the vortex. Most of the 
time vortices only appear with topological charges of ±1. The two signs represent two opposite helicities, which 
indicate whether the phase increases in a clockwise or counterclockwise direction around the phase singularity. 

The rate of increase of the phase around a vortex is not always a constant function of the azimuthal angle. 
For a canonical vortex the rate of increase is constant and the vortex has an isotropic (rotationally symmetric) 
appearance. One can regard a noncanonical vortex as the anisotropically scaled (or stretched) version of a 
canonical vortex. Together with the direction or orientation of the scaling, the amount of anisotropy represents 
the morphology of the optical vortex. 

One can create an optical vortex in a beam with a computer generated hologram. 3 This vortex then propagates 
along the beam to the point where it can be applied to rotate some small objects. 4 Unfortunately, the vortex 
does not always remain the same as it propagates. If there are more than one vortex in the beam they interact 
with each other and can even annihilate each other. 5 Moreover, in the late 90's a group in spain 6 performed 
an interesting experiment in which they demonstrated that the topological charge of an optical vortex can be 
inverted after the vortex bearing beam passed through a cylindrical lens. Hence, under certain circumstances, 



even if there is only one vortex in the beam, it can change its shape (morphology) and even its topological 
charge. So the topological charge in a beam is not always conserved during propagation. What happens is that 
two 'vortices' are brought in from infinity by means of an edge phase dislocation to cancel the vortex in the beam 
and replace it with a vortex of opposite topological charge. This happens because of the astigmatism that the 
beam picked up from the cylindrical lens. The topological charge is still locally conserved in that the net inward 
(or outward) flow of topological charge through a closed surface is always zero. 

One of the consequences of this observation is that there is a set of special planes in such an astigmatic beam 
where the topological charge is inverted. 7 One of the questions one can ask is where these planes are located 
and what these locations depend on. Here we report on the results of this investigation. It is well known that 
a single noncanonical vortex that propagates in a rotationally symmetric, stigmatic Gaussian beam retains its 
morphology for as far as it propagates. 8 When the rotational symmetry is broken or the beam picked up some 
astigmatism this situation is no longer true - the morphology of the vortex evolves during propagation. If the 
Gaussian beam is asymmetrical (elliptical) but not astigmatic and the vortex is canonical in the waist of the 
Gaussian beam, then, although the vortex becomes noncanonical, it does not change its topological charge. If it 
is noncanonical in the waist, there can be a region over which the topological charge of the vortex is inverted. A 
more detail analysis of these conditions is provided below. If the beam is astigmatic but symmetric in the sense 
that the Rayleigh distances associated with the two transverse directions are equal, then the vortex can invert 
its topological charge irrespective of its initial morphology. Hence it is both the morphology of the vortex and 
the astigmatism of the beam that determines whether topological charge inversion occurs. 

The rest of the paper is organized as follows. We define our notation for the beam and the morphology of the 
vortex in Section 2. The general expressions for the trajectory of a single vortex with an arbitrary morphology 
propagating in a general asymmetrical and/or astigmatic Gaussian beam are provided in Section 3. The zeros 
of the denominator of these trajectories provide the locations for the flip planes. In Section 4 these flip plane 
locations are investigated for various situations. We provide expressions for the anisotropy of the vortex as a 
function of the propagation distance in Section 5 to show that the anisotropy vanishes and changes sign at the 
locations of the flip planes, which confirms that topological charges inversion occurs at the flip plane locations. 
Conclusions are provided in Section 6. 

2. NOTATION 

A vortex bearing Gaussian beam can be expressed as a product of a complex valued polynomial function (i.e. the 
prefactor) and a Gaussian function. We'll define the propagation direction as the z-direction and the transverse 
directions as the x- and y-directions. Here, we'll use normalized coordinates u = x/ujq, v — y/uJo and t = z/p, 
where wo and p represent the waist size and the associated Rayleigh range of the Gaussian beam, respectively. 

An astigmatic Gaussian beam has two different waists (foci) for two orthogonal transverse directions, which we 
define as the x- and ^/-directions, respectively. We define a central plane as the plane that is exactly halfway 
between the two focal planes of the astigmatic beam. For a stigmatic beam the central plane coincides with 
the waist of the Gaussian beam. One can then quantify the amount of astigmatism by a parameter r, which is 
defined as the distance from the central plane to one of the focal planes, normalized with respect to the nominal 
Rayleigh range (defined below). Note that r is a dimensionless parameters that is given as the ratio of the two 
most significant scale parameters in the beam. 

A Gaussian beam can be asymmetric in the sense that the size of the beam waist may be different for the x- and y- 
directions, respectively. In this case we define, for the purpose of normalization, a nominal waist size luq = ^Ju) x uj y , 
where u> x and u> y are the waist sizes in the x- and y-directions, respectively. The associated nominal Rayleigh 
range is then p = ^Jp x p y . We quantify the asymmetry of the beam by a parameter 7 = lu x /lu v = p x /p y - For 
definitencss we'll assume that 7 > 1 u x > ui y . For 7=1 the beam is called symmetric (which is not to be 
confused with rotationally symmetric beams where we also have r = 0). 



The prefactor for a Gaussian beam with a single on-axis noncanonical vortex is given by 9 

V(u,v) = -j=[S(u + iv) + C(u-iv)], (1) 

where the morphology of the vortex is given in terms of 10 

£ = cos(V>/2) exp(-i<f>/2) (2) 
C = sin(^/2)exp(^/2), (3) 

with tp and 4> being the morphology angles, representing the anisotropy and orientation, respectively. The sign 
of cosV> represents the topological charge of the vortex. For an off-axis vortex the expression in Eq. (1) can be 
shifted to an arbitrary location (uo,Vo). 

Hence, in its central plane a general asymmetric astigmatic Gaussian beam with a single off-axis noncanonical 
vortex can be expressed by 



9{x,y) = -j= [(£ + C)(«-«o) +i{£-0(v -v Q )]exp ^- 
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The expression at other points along the direction of propagation can be obtained with the Fresnel transform 11 
of this expression. 

3. TRAJECTORIES 

The prefactor changes as a function of t. The zeros of the prefactor give the locations of the vortex as a function 
of the propagation distance - the vortex trajectories. For a general asymmetric astigmatic Gaussian beam with 
a single off-axis noncanonical vortex, the vortex trajectory is given by 

( ] = [ 7 2 + (r - tf] {7 [(Ct 7 - A)u + (1 - B)v ] t + C(tV + l)u } 

U[ ' L4r 7 (l + 7 2 ) + C(l - r 2 ) 7 2 ] t 2 - [A 7 (7 2 - 1)(1 - r 2 ) - Cr( 7 4 - 1)] t + C(t 2 ^ 2 + l)(r 2 + 7 2 ) ' U 

( ] - [1 + (t + r) 2 7 2 ] { [(1 + B)>yu + (Ct - A^vp] t - C(t 2 + ^vp } 

V[ ' [Aty(1 + 7 2 ) + C(l - r 2 ) 7 2 ] t 2 - [A 7 ( 7 2 - 1)(1 - r 2 ) - Cr( 7 4 - 1)] t + C(T 2 ~f 2 + l)(r 2 + 7 2 ) [ ' 
where 

A = sin ipo sin <f) B = sin ip cos <po C = cos V'o, (7) 
with ip and $0 being the morphology angles in the central plane. 

4. FLIP PLANE LOCATIONS 

The values of t for which the denominator of the trajectory functions, given in Eqs. (5-7), become zero represent 
points along the propagation direction where the vortices move to transverse infinity. As we'll show below the 
topological charges are inverted in the planes at these locations. These planes are therefore called flip planes. 
Since the denominator is quadratic in t there arc in general two such planes in a beam. The existence of these 
planes in astigmatic Gaussian beams is in itself an interesing fact, because they do not in general coincide with 
any of the other special planes in the beam, such as the focal planes. 

Note that the denominator in Eqs. (5) and (6) is independent of the coordinates of the vortex in the central 
plane, (u 07 v ). It only depends on the astigmatism parameter r, the symmetry parameter 7 and the initial 
morphology of the vortex. 



We now consider the locations of the flip planes under two special circumstances: for elliptical stigmatic beams 
(t = 0) and for symmetric astigmatic beams (7 = 1). 



4.1. Elliptical stigmatic beams 

First we consider the location of the flip planes in an elliptical stigmatic beam and thus set r = in the 
denominator of the trajectory functions. The locations of the flip planes are then given by 



where 



and 




Figure 1. The locations of the flip planes in a stigmatic beam as a function of K. 

The locations of the flip planes as a function of K are shown in Fig. 1. One can see that the flip planes are 
always located on the same side of the central plane. In the extreme case (K — > oo) when the vortex is an edge 
phase dislocation there is only one flip plane located in the central plane. 

We see that there are only flip planes when 

2 7 

| tan ip sin <f> \ > -5 -, (11) 

— I 

which implies that a certain amount of anisotropy is required for the vortex in the initial plane to produce 
flip planes. A canonical vortex in the waist cannot undergo topological charge inversion in a stigmatic beam. 
Moreover, if the beam is symmetric (7 = 1) in addition to being stigmatic (r = 0) — in other words, the beam is 
rotationally symmetric - then a finite amount of anisotropy cannot satisfy the requirement of Eq. (11). 



4.2. Symmetric astigmatic beams 

Now we consider the location of the flip planes in a symmetric astigmatic beam for which we set 7 = f in the 
denominator of the trajectory functions. In this case the locations of the flip planes are given by 



ttp=± 



1 



<t« - 2Tt - 1 

with T given in Eq. (10). We see that the flip planes only exist when 



t>T+ VT2 + 1. 



(12) 



(13) 



Note that the two flip planes are now located at equal distances from the central plane and they lie on opposite 
sides of the central plane. 



The minimum flip plane distance for each value of r occurs when 



T = tan "00 sin 4>q 



r(r 2 -3) 
3r 2 -l • 



and this minimum distance is given by 
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(14) 



(15) 



This implies that the flip plane moves all the way to the central plane (t = 0) when r = 1/V3, for which T — > 00. 
In other words, the vortex in the central plane becomes an edge dislocation. 
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Figure 2. The flip plane distance in a symmetric beam (7 = 1) with a canonical vortex in the central plane (T = 0), as 
a function of r. 

The flip plane distance is shown in Fig. 2 as a function of r for the canonical case (T = 0). Here flip planes exist 
when r > 1. The flip plane distance starts at infinity for r = 1, reaches a minimum at r — \/3 and rises again 
to approach infinity for r — > 00. The minimum flip plane distance at r = v3 is t = 2y/2. 



5. ANISOTROPY 



The significance of the existence of the flip planes lies in the fact that the vortices change their topological charge 
in those planes. One can show this by computing the anisotropy of the vortex as a function of the propagation 



distance. The anisotropy of a vortex in a complex valued function / is given by 5 

\d-f\ 2 -\d + f\ 2 



cosi/j 



evaluated at the location of the vortex, where 



d ± = — (— ±i— 

y/2 \dx dy 



(16) 



(17) 



For an asymmetric astigmatic Gaussian beam with a vortex of arbitrary morphology in the central plane, the 
anisotropy of the vortex as a function of the propagation distance is given by 



2 { [An(l + 7 2 ) + C(l - r 2 ) 7 2 ] t 2 - [Aj(j 2 - 1)(1 - r 2 ) - Cr( 7 4 - 1)] t + C(t 2 j 2 + l)(r 2 + 7 2 )} 



(18) 



[B( 7 4 - 1) + 7 4 + 2 7 2 r 2 + 1] t 2 + 2r [B( 7 4 + 2 7 2 r 2 + 1) + 7 4 - 1] t + 2(1 + 7 2 r 2 )( 7 2 + r 2 ) 

In the central plane the anisotropy is given by cosV'o- From the fact that the numerator of the anisotropy 
function in Eq. (18) is proportional to the denominator of the trajectories, provided in Eqs. (5) and (6), one can 
conclude that the anisotropy becomes zero. The topological charge is given by the sign of the anisotropy. One 
can therefore see that the topological charge changes sign at the flip planes. 

Next the anisotropy is considered for elliptical stigmatic beams and for symmetric astigmatic beams. 
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Figure 3. The evolution of the anisotropy (cosip(t)) for a stigmatic asymmetric beam with j — 4 and a vortex morphology 
given by ipo = 4>o — 7r /4- 



5.1. Elliptical stigmatic beams 

The anisotropy for an elliptical stigmatic beam is given by 

2 [cosi/'o7 2 i 2 — sin^o sin(/)o 7 ( 7 2 — l)t + cos^o 7 2 l 



cos ip(t) = 



[sinV-o cos 0o (7 4 - 1) + 7 4 + 1] t 2 + 2 7 2 
If the vortex is canonical in the central plane (tpQ — 0, 7r), the anisotropy becomes 

±2 7 2 (t 2 + 1) 



cos tp(t) 



(7 4 + l)i 2 + 2 7 2 



(19) 



(20) 



Apart from the sign in front of it, the right-hand expression in Eq. (20) is positive, which means that there are 
no flip planes. It therefore requires a non-canonical vortex in the central plane to produce flip planes. To find 
the minimum required morphology, one can consider where the numerator of Eq. (19) becomes zero and find the 
same requirement given in Eq. (11). 

In Fig. 3 we show the evolution of the anisotropy for a beam with r = and 7 = 4 and a vortex with morphology 
given by tpo — (f>o — it /A. One can see that there is a range of propagation distances on one side of the central 
plane for which the anisotropy is negative, hence, where the topological charge of the vortex is inverted. We also 
see that the maximum anisotropy is achieved in a plane different from the central plane. 



5.2. Symmetric astigmatic beams 

The anisotropy for a symmetric astigmatic beam is given by 

2rsmVosin0 o i 2 + cosf/'o [(1 - T 2 )t 2 + (1 + r 2 ) 2 



cosif)(i) 



(1 + t 2 + t 2 + 2rsmip cos</> i)(l + t 2 ) 



(21) 



For ?/>o ^ the vortex is not canonical in the central plane. However, it becomes canonical at another point 
along the propagation direction if 

COS -0Q ± 1 

T = ; — ; — ; — r~, (22) 

sm yjQ sin ipo 

where the sign corresponds to the sign of the topological charge of the vortex in the central plane. The location 
where the vortex becomes canonical is 

< + = — tan c 



1 =F cos 0o 



(1 ± cos ipo ) cos </>o sin 0o ' 

where the upper (lower) sign represent a canonical vortex with a positive (negative) topological charge. 



(23) 



o 

o 
.52 
'E 
< 




-10 



-6 -4 -2 2 4 6 8 

Propagation distance [normalized units] 



10 



Figure 4. Anisotropy of a canonical vortex in a symmetric astigmatic beam with r = 2, as a function of propagation 
distance t. 

The anisotropy is shown in Fig. 4 as a function of t for the canonical case with r = 2 and 7=1. The function 
is symmetric around the y-axis because it only depends on t 2 . At the central plane the anisotropy is 1, which is 



the value with which it has been launched. This value decreases and reaches zero at the flip planes. For t — > oo 
the anisotropy approaches the value —3/5. 

6. CONCLUSIONS 

We investigated the properties of the planes in astigmatic and/or elliptical beams where the topological charge 
of optical vortices are inverted. On the one hand, these flip planes are obtained as the zeros of the denominator 
of the trajectory equations. On the other hand, they are given as the zeros in the anisotropy of the vortex as a 
function of the propagation distance. Two special cases are considered: elliptical stigmatic beams and symmetric 
astigmatic beams. For the first case flip planes only exist if the vortex is sufficiently anisotropic in the waist 
(central plane) of the beam. In the second case flip planes exists even for vortices that are canonical in the 
central plane, provided that the beam is sufficiently astigmatic. 

The behavior of vortices and in particular the location of the flip planes in astigmatic and/or elliptical beams 
can thus be predicted. This provides important information for the design of optical systems in the application 
of optical vortices. 
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